Laser-Generated Ultrashort Multi-Megagauss Magnetic Pulses in Plasmas 
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We demonstrate ultrashort (6 ps), multi-Megagauss (27 MG) magnetic pulses generated upon 
interaction of an intense laser pulse (10 16 Wcm -2 , 100 fs) with a solid target. The temporal 
evolution of these giant fields generated near the high density critical layer is obtained with the 
highest resolution reported so far. Particle-in-cell simulations and phenomenological modeling is 
used to explain the results. The first direct observations of anomalously rapid damping of plasma 
shielding currents produced in response to the hot electron currents penetrating the bulk plasma 
are presented. 
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The largest magnetic fields available terrestrially (~ 
10 8 Gauss) are generated by explosive ionization of a solid 
target with an intense ultrashort laser pulse yj. Since 
the first observation of such magnetic fields, their origin, 
magnitudes and other qualitative features have attracted 
considerable attention [ 2) . Recently, sub-picosecond laser 
produced solid plasmas have provided a new experimen- 
tal facet to these studies. Magnetic fields up to gigagauss 
magnitudes have been predicted in overdense region of 
solid target j|. Little, however, is known about the tem- 
poral evolution of these huge magnetic fields generated 
around critical layer. These fields play a crucial role in 
electron transport jj] and are, therefore, important for 
potential applications in hybrid confinement Q and fast 
ignition g] schemes of laser fusion. 

In this Letter, we present first experimental measure- 
ments of temporal evolution of Megagauss magnetic fields 
generated at the critical layer, on femtosecond timescales. 
The field generation and decay mechanisms are identified 
and role of resonance absorption is examined. The initial 
buildup of magnetic field due to direct laser radiation ef- 
fects is calculated using LPIC++ code and the results 
are found to be in good agreement with the experiment. 
The field evolution is explained to be due to currents 
generated by fast electrons ]8| and plasma return cur- 
rents damped by turbulence induced resistivity ||. The 
first direct observation of anomalously rapid damping of 
return currents, which may have important consequences 
for laser fusion, is reported. We demonstrate ultrashort, 
Megagauss magnetic pulses, with 6 picosecond (FWHM) 
duration and a peak magnitude of 27 Megagauss gener- 
ated by a p-polarized laser pulse. Possible applications 
of such high-field ultrashort magnetic pulses in various 
areas are discussed. 

Laser pulses (806 nm, 100 fs) are derived from a cus- 
tom built chirped pulse amplification Titanium:Sapphire 
laser, described in detail elsewhere |l(|. The linearly po- 
larized pump pulse incident at 55° is focused on a solid 
target with typical intensities of 10 16 W cm -2 (Fig. 1). 
A small part of the laser is used to generate a linearly 



polarized, second harmonic (403nm) probe pulse using a 
thin (0.5mm) BBO crystal. The probe pulse incident at 
50°, is focused to an intensity of 5 x 10 12 W cm~ 2 and 
is spatially overlapped with the pump spot. The probe 
penetrates beyond the critical density(n c ) for pump (i.e. 
upto 1.65n 806 or 0.4n* 03 as shown in inset of Fig.l). The 
probe pulse delay with respect to pump is varied using a 
motorized translation stage. 

The novel features of this experiment are: (a) we are 
able to examine the generated magnetic field on either 
side of the critical layer by using the second harmonic of 
the laser as a non-tangential probe, (b) we identify zero- 
delay time precisely and obtain high-resolution data for 
temporal evolution, and (c) we examine, using oblique 
pump incidence, the role of Resonance Absorption (RA). 
To our knowledge, only one laboratory has hitherto re- 
ported |llj] pump-probe magnetic fields measurements 
on picosecond time scales (using 1 ps duration laser 
pulse at normal incidence). However, these measure- 
ments were confined to an under-dense plasma region 
(n e = 4 x 10 19 cm~ 3 ) and the diagnostics employed lim- 
ited the temporal resolution to ~ 3 ps. We emphasize 
our probing technique, because, magnetic field genera- 
tion primarily occurs near the critical surface (electron 
density ~ 10 21 cto~ 3 ), the region of maximum laser ab- 
sorption ||. Moreover, it is the magnetic field in the 
over-dense region that determines hot electron transport 
into the bulk, which is crucial for fusion related issues. 

We use standard Cotton-Mouton polarimetry |fl2| , 
which involves measurement of magnetically induced el- 
lipticity in the polarization state of the probe pulse. The 
contributions to induced ellipticity from plasma propa- 
gation and refraction effects |ll| are calculated and are 
found to be insignificant |Q. The targets used in these 
experiments are optically polished aluminum discs. The 
target is housed inside a vacuum chamber at 10 -3 Torr 
and is translated after each laser shot, exposing a fresh 
surface for interaction. The pump and probe focussing is 
monitored using hard x-ray emission as diagnostic. The 
spot radii for pump and probe are 15//m and 10/^m re- 
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FIG. 1: The experimental setup, Target (T), Lens (L), Beam 
Splitter (BS), Photo Diode (PD), Half Wave Plate (HWP). 
Inset shows diagrammatic representation of pump and probe 
laser paths in the plasma. The critical densities for pump and 
probe are given as nj? 06nm and nf 3nm respectively 



FIG. 2: Magnetic field pulse profile for p- and s-polarized 
pump (open squares and triangles respectively) laser with in- 
tensity of 1.1 x 10 16 W cm" 2 . Solid line shows the fit ob- 
tained in case of p-polarized pump using our phenomenologi- 
cal model. The inset shows the reflectivity and induced ellip- 
ticity of probe as function of delay time. 



spectively. The reflected probe beam is collected using 
f/20 lens. The amplitude and polarization of the specu- 
larly reflected probe is simultaneously monitored by split- 
ting it into two arms, one measuring the reflectivity with 
a photodiode (PD2) and the other measuring the polar- 
ization state using an analyzer (extinction ratio 10~ 5 ) in 
front of another identical photodiode (PD3), as shown in 
Fig. 1. The photodiode PD1 measures the shot to shot 
laser fluctuations. All three PD signals are simultane- 
ously recorded for each laser shot along with the delay 
stage position. The ellipticity is determined at various 
fixed delay positions by plotting PD3 signal with respect 
to analyzer angle. Higher temporal resolution elliptic- 
ity data is obtained by continuously varying the delay 
in steps of lfj,m at different fixed analyzer positions and 
computing the ratio PD3/PD2 for each delay. The ratio 
(PD3/PD2) is used above so as to account for plasma 
reflectivity variation as function of time delay. The tem- 
poral behavior of the magnetic field from -2 to +10 ps 
time delay is monitored. The constant background noise 
due to the second harmonic radiation from the plasma 
generated by the pump pulse is negligible (< 0.1%) as 
compared to the reflected probe. 

Figure 2 presents the temporal evolution of the mag- 
netic field. The inset shows reflectivity and induced ellip- 
ticity of the probe for a p-polarized pump. The sharp re- 
flectivity dip is used to independently establish the start 
of the magnetic pulse. The magnetic field is derived 
from induced ellipticity (/3) for our experimental condi- 
tions using H f3(t) = 3.32 x 10" 26 / n e (l,t)B 2 (l,t)dl, 



where 'n e ' is the electron density in cm -3 , 'B' is mag- 
netic field in MG, and T is the path length in /nm. The 
magnetic field is deduced assuming a spatially uniform 
B over a linear density gradient n e (x) = 0.4n* 03 (a;/£), 
where L is the plasma slab length, as shown in the inset 
of figl . The factor 0.4n^ 03 corresponds to the turning 
point density for the 403nm probe. We integrate over 
the trajectory in the plasma, dl = yl + y' 2 dx, where 
y' = sin9o/ y/e(x) — sin 2 0Q , e(x) being the dielectric 
function. The plasma expansion velocity is estimated 
to be 5 x 10 7 cm/s from Doppler shift measurements of 
the reflected probe. With L{t) = 1 + v exp t , the mag- 
netic field as a function of time delay is obtained as 
B{t) = 80v7?(*)/(l + 0.5f) MG . From the results shown 
in Fig. 2, the magnetic field pulse generated p-polarized 
pump has a peak value of 27 MG and duration (FWHM) 
of 6ps. In comparison s-polarized pump results in a peak 
value of 14 MG. This confirms the importance of RA, 
induced by p-polarization, in magnetic field generation. 
The magnetic field in case of s-polarized pump can be 
explained to be arising from parametric instabilities near 
critical density, which develop at a slower rate as seen in 
figure 2. Further, any realistic laser focusing always in- 
volves some RA contribution even for s-polarization, due 
to critical surface rippling and geometrical effects || [l5) • 

The mechanism of quasi-static magnetic field genera- 
tion in ultrashort laser -plasma experiments can be un- 
derstood by invoking the following magnetic field evolu- 
tion equation 
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Equation (1) is derived by taking the curl of the equa- 
tion of motion of background plasma electrons that carry 
the plasma shielding currents, j p — V x B — jhot ■ The 
first term in equation (1) is the standard electron mag- 
netohydrodynamic (EMHD) Jl6[ source due to the Hall 
effect, the second term is the thermoelectric source and 
the third term is the source due to hot electrons (gen- 
erated near the pump critical surface by the RA of the 
laser pulse); the last term gives the magnetic field decay 
due to resistive damping of the plasma shielding currents 
(ct -1 being the background plasma resistivity). Electron 
inertia and magnetic field convection effects are assumed 
to be negligible in equation (1). 

We estimate the first term (important during the first 
200 fs when the pump is on) by taking a cycle-averaged 
product of high frequency current and magnetic fields 
p7[,V x (^(j h x B h )/n e e^, inside the plasma. We obtain 
the magnitude of this term by carrying out a particle- 
in-cell simulation using the laser-plasma interaction code 
LPIC++. In the simulation, a p-polarized light pulse 
with sin 2 envelope and 100 fs FWHM laser pulse is in- 
cident at an angle of 55° on a linear density ramp. Fig. 
3(a) shows numerically obtained spatial and temporal 
profiles of the quasi-steady state magnetic field. This 
confirms the earlier analytical results of spatial profiles 
of laser generated magnetic fields ]Is[ ]. At the same time, 
the values of B and dB/dt obtained from simulation 
compare well with our experimental results (Fig. 3b). 
The numerically obtained maximum value of dB/dt is 
22 MG/ps at 120 fs, which is close to the experimentally 
deduced dB/dt. The B value reached in 150 fs is 1.2 MG, 
beyond which the simulation results saturate. This is ex- 
pected, because to model this regime we have considered 
only one term, which is valid strictly for the duration of 
pulse. 

The second term in equation (1) is estimated assum- 
ing a temperature gradient of 100 eV over a transverse 
scale of 15/Ltm (pump spot radius) and density gradient of 
10 21 cm -3 over 1/^m in the normal direction. This gives 
dB/dt ~ 0.05 MG/ps, which is much smaller than the 
experimental results. Hence, the thermoelectric source 
term is neglected in subsequent analysis. 

The remaining terms (magnetic diffusion and hot elec- 
tron source ^, |l{| ) in the magnetic field evolution equa- 
tion govern the evolution of B after the pump pulse is 
removed. We model this regime by a phenomenological 
0-dimensional evolution equation 
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FIG. 3: (a) The temporal and spatial profile of magnetic field 
obtained using the LPIC++ code (x = x c is the critical layer), 
(b) Comparison of B and dB/dt obtained using LPIC++ sim- 
ulation with experimental results. 



the fast electron currents and B/t is a 0-d represen- 
tation of the magnetic diffusion term. Taking S(t) = 
Soexp(—t/to), we get 



B(t) 



S'o 



1/t - 1/t, 



■[exp(-t/t ) - exp(-t/r)], (3) 



where S(t) is the source term. S(t) is mainly due to 



The first term in equation (3) denotes the natural decay 
of the hot electron source produced by RA mechanism. 
The second term describes the resistive decay of the fields 
generated by the plasma return currents. As shown in 
Fig. 2, this expression gives an excellent fit to our ex- 
perimental data for p -polarized pump with S'o = 53.7 
MG/ps, t = 0.7 ps and r = 5.6 ps. 

To get an insight into these numerical values, we es- 
timate conductivity a from the magnetic diffusion term. 
Using r ~ (4tt(t/c 2 )(Ax) 2 , for the best fit value r = 5.6 
ps and Ax = 15/im , we get a — 1.8 x 10 14 sec -1 . How- 
ever the <J solid Or <J classical observed [§0) at T e = 100 
eV has value of rs 4.5 x 10 15 sec -1 , which is an order 
of magnitude greater than that obtained from magnetic 
field decay! This clearly brings out the importance of 
turbulence induced anomalous resistivity effects in the 
damping of shielding plasma currents. An upper esti- 
mate of anomalous resistivity due to turbulence of elec- 
trostatic waves, obtained by taking the effective collision 
frequency v e ff ~ f^>p ( where / is a fraction of order 
unity ), is in reasonable agreement with our estimate of 
a. We note that the hot electron source term is effec- 
tive for about 0.7 ps, which is longer than the laser pulse 
(FWHM=100fs). Thus we argue that the electrostatic 
plasma waves generated by the RA mechanism (typically 
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with E 2 /4irn e T e ^> 1 ) during the laser pulse will con- 
tinue to slowly damp (or convect away) and accelerate 
electrons for a few hundred femtoseconds after the laser 
pulse. To get an estimate of So, which is the magni- 
tude of the source term | (c/a)(V x jhot) \ , we evaluate 
jhot ~ f a eI/T hot by taking a conversion fraction (/ a ) of 
incident energy into hot electrons as 0.3 and That ~ 20 
keV (estimated by using the well known scaling laws for 
resonance absorption This yields hot electron cur- 

rent density jhot ~ 4.5 x 10 20 statampere/cm 2 . Using 
anomalous conductivity obtained above we get So ~ 61 
MG/ps, which is again in close agreement with our phe- 
nomenological fit. 

In conclusion, we have measured and characterized pi- 
cosecond megagauss magnetic pulses generated by the in- 
teraction of ultrashort laser pulse with a solid. Our mea- 
surements extend to overdense region of the target and 
hence are of relevance to electron transport and fusion 
related issues. The experimentally observed rise times 
and magnitude of magnetic fields closely follow theoret- 
ical estimates and simulations. We also observe for the 
first time, anomalously rapid damping of return plasma 
shielding currents produced in response to the hot elec- 
tron currents penetrating the bulk plasma; this is a topic 
of great significance to the fast ignition scheme. Such 
ultrashort, localized magnetic fields are useful for in- 
vestigating magnetic precession and reversal dynamics, 
which is vital for developing next generation ultrafast 
switching and storage devices p2| |. Further, the genera- 
tion and characterization of these giant magnetic fields 
offers a unique opportunity for accessing extreme stellar 
conditions and testing astrophysical theories in the lab- 
oratory j2^]. The diagnostics in chemical and biological 
sciences like Magnetic Circular Dichroism and Magnetic 
Resonance may also benefit in situations where high fields 
are essential |24| . The intensities used in our experiments 
are easily realizable with modern kilohertz repetition rate 
femtosecond lasers, and we foresee exciting applications 
for these magnetic pulses. 
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